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Abstract 

The presence of several clone solutions in the simultaneous measurement of 
($13> 5) has been widely discussed in literature. In this letter we write the analytical 
formula? of the clones location in the (#13, <5) plane as a function of the physical input 
pair (0i3, 8). We show how the clones move with changing #13. The "clone flow" can 
be significantly different if computed (naively) from the oscillation probabilities or 
(exactly) from the probabilities integrated over the neutrino flux and cross-section. 

Using our complete computation we compare the clone flow of a set of possi- 
ble future neutrino experiments: the CERN SuperBeam, BetaBeam and Neutrino 
Factory proposals. We show that the combination of these specific BetaBeam and 
SuperBeam does not help in solving the degeneracies. On the contrary, the combi- 
nation of one of them with the Neutrino Factory Golden and Silver channel can be 
used, from a theoretical point of view, to solve completely the eightfold degeneracy. 
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1 Introduction 



The atmospheric and solar sector of the PMNS leptonic mixing matrix pQ have been 
measured with quite good resolution by SK [2|, SNO [3] and KamLand jlj. These ex- 
periments measure two angles, 9 12 and #23, and two mass differences, Am 2 2 and Am^ 
(for the explicit form of the PMNS matrix and the adopted conventions, see for example 
jSJ). The present bound on 6>i 3 , sin 2 #13 < 0.02, is extracted from the negative results of 
CHOOZ 6j and from three-family analysis of atmospheric and solar data. The PMNS 
phase 8 is totally unbounded as no experiment is sensitive to the leptonic CP violation. 
The main goal of next neutrino experiments will be to measure these two still unknown 
parameters of the leptonic mixing matrix. In this paper we concentrate on analyzing, 
form a theoretical point of view, a particular problem that arise when trying to measure 
simultaneously (0i 3 , S). 

In [7] it has been noticed that the appearance probability P a p(9i3,S) obtained for 
neutrinos at a fixed energy and baseline with input parameter (#13, 5) has no unique 
solution. Indeed, the equation: 



has a continue number of solutions. The locus of (6*13, 5) satisfying this equation is called 
"equiprobability curve" . Considering the equiprobability curves for neutrinos and antineu- 
trinos with the same energy (and the same input parameters), the system of equations 



(where ± refers to neutrinos and antineutrinos) has two intersections: the input pair 
(6*13, 5) and a second, energy dependent, point. This second intersection introduces an 
ambiguity in the measurement of the physical values of 9^ and 5: the so-called intrinsic 
clone. As it was made clear in [HI IH] 5 two other sources of ambiguities are present: 

• We only know the absolute value [2|, |Am| 3 | G [1.9,3.5] x 10~ 3 eV 2 of the atmo- 
spheric mass difference, not its sign. This ambiguity originates from the fact that 
atmospheric neutrino experiments measure only the leading oscillation — > v T that 
depends quadratically on Am| 3 ; 

• We only know the allowed departure of #23 from maximal mixing, sin 2 19^3 > 0.9, but 
not if #23 is actually smaller or greater than 45°. This ambiguity originates in that 
neutrino experiments are looking for v e disappearance or — > v T oscillation, 
where #23 appears only through sin 2 2# 2 3- 

As a consequence, future experiments will have as ultimate goal the measure of the 
two continuous variables #13 and 5 plus the two discrete variables: 



P a p(9i3,6)=P a p(9i3,5) 




P±(9 13 ,6)=P^(9 13 ,5) 



(2) 



Satm = sign[Aml 3 ] , 
Sod = sign[tan(29 2 3)] ■ 



(3) 
(4) 
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These two discrete variables assume the values ±1, depending on the physical assignments 
of the Am| 3 sign (s atm = 1 for m\ > m\ and s atm = — 1 for m\ < m|) and of the #23- 
octant (s oct = 1 for 2 3 < tt/4 and s oct = —1 for 9 23 > 7r/4). As a consequence, taking 
into account all the present ignorance on the neutrino mixing and masses, eq. (J2J) must 
be rewritten, more precisely, as: 

Pap(.@13i ^] Satm-, Soct) — P a p(@137 8] Satm = S~atm] S oc t = Soct) , (5) 

where s atm and s oct have been included as input parameters in addition to 6> 13 and 6. 
In eq. (JSJ) we have implicitly assumed to know the right sign and the right octant for 
the atmospheric mass difference and angle. As these quantities are unknown (and pre- 
sumably they will remain so at the time of the next neutrino facilities), all the following 
equiprobabilities systems of equations should be considered as well: 

Papi^^'it 8] Satm-, Soct) = Pap{@lZi S a t m = —Satm] S oc t = Soct) , (6) 
^a^(^13) 8] S a tm, Soct) = P a p{^^i ^'i s atm = Satm] Soct = ~Soct) i (7) 
Pap{@13i 8', S atrn , S oct ) = Pafi(9i 3 , 6] S atm = —Satm] S oc t = ~Soct) ■ (8) 

These new sets of equiprobability systems arise when we equate the true probability 
(left hand side) with the probabilities obtained including one of the three possible wrong 
guesses on s atm and s oct (right hand side). 

Solving the four systems of eq. ©-© will result in obtaining the true solution plus the 
appearance of additional clones to form an eightfold-degeneracy 0. These eight solutions 
are respectively: 

• the true solution and its intrinsic clone, obtained solving the system of eq. (0); 

• the Am^-sign clones (hereafter called sign clones) of the true and intrinsic solution, 
obtained solving the system of eq. (jUJ); 

• the #23-octant clones (hereafter called octant clones) of the true and intrinsic solu- 
tion, obtained solving the system of eq. (J7J); 

• the Am^ m -sign 6* 2 3-octant clones (hereafter called mixed clones) of the true and 
intrinsic solution, obtained solving the system of eq. (jSJ). 

In the following section we illustrate the technique used for calculating the location 
of the clone solutions. We describe explicitly the formulae for the v e — ► oscillation 
probability, due to the great importance it has in literature. In fact this "golden" channel 
(together with its CP conjugate channel) has been shown to be very promising to simulta- 
neously study the 9% 3 angle and the leptonic CP-violation at the Neutrino Factory [TU1 E] 
and in the BetaBeam based facilities ^2]- Anyway most of the considerations done apply 
as well to the time reversal transitions — > u e (and the CP conjugate one) that has been 
considered in the context of SuperBeam facilities like the SPL CERN project [12]. The 
same strategy will be applied in Sect. (j3.3|) to the study of the v e — > v T transition (the 
so-called "silver channel" ) at a Neutrino Factory. As pointed out in ^3 this channel 
can become very useful in reducing the number of clone solutions. 
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2 Probability vs Number of Events 



Following eq. (1) of [7j, the appearance probability for the v e — > {p e — > z^J transition, 
expanded at second order in perturbation theory in 9 13 , A 12 /A 2 3, A 12 /A and A 12 L (see 
also [HUE]), reads: 

P±(0 13 , 5) = X ± sin 2 (2# 13 ) + Y± cos(e 13 ) sin(2# 13 ) cos (±5 - ^) + Z , (9) 

where ± refers to neutrinos and antineutrinos, respectively. The coefficients X^, Y and 
Z are defined as follows: 




with A = \f2Gpn e (expressed in eV 2 /GeV) and B T = \A =F A 23 |. Finally, A 2 3 = 
Am 2 l 3 /2E L/ and A 12 = Am\ 2 /2E U . Instead of using the coefficients in eqs. (jlOH12j) . it 
is convenient [TH] to introduce the following (#i 3 , 5)-independent terms: 

Of = X\ 0± = y±cos(^), Ot = Y ± mn(^) f A = Z. (13) 

Notice, however, that transition probabilities, like for example P^, are not the ex- 
perimetally "measured" quantities. Experimental results are given in terms of the number 
of charged leptons observed in a specific detector. For the Neutrino Factory Golden chan- 
nel, for example, one measures the number of interacting muons with charge opposite to 
those circulating in the storage ring. In a SuperBeam-driven detector the signal is, on the 
other hand, represented by electrons. If the considered detector can measure the lepton 
and hadron final energies, Ei, Eh with a good resolution, the events are then grouped in 
bins of energy interval AEi. This is indeed the case for the proposed Magnetized Iron 
Calorimeter (MIC) 19j at the Neutrino Factory, where the number of muons in the i-th 
energy bin for the input pair (# 13 , 5) and for a parent muon energy E^ is given by: 

N^(9 13 ,5) - | — ® P e ,(E v ,9 13 , 5) ® — 1^ (14) 

where <S> stands for a convolution integral [T3|. The same happens for the Emulsion Cloud 
Chamber (ECC) (see Sect. 13 .31 for details). On the contrary, whenever the considered 
detector energy resolution is not sufficient to group events in bins of energy, one single 
sample will be considered. This is for example the case for the Water Cherenkov detector 
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located at the Frejus as the target of a SuperBeam (SB) or Beta-Beam (BB) experiments, 
as we can read from 



Using the Oj coefficients in eq. (|T3j) we define the following integrals: 

Ei+AE/j 

3 (E u ) ® ][r - r '" ) . 1 



dE„ 



dE tJ 



(15) 



where the Oj coefficients are convoluted over the vN cross-section and the differential 
neutrino flux. All the experimental details, such as the flux shape, the cross-section, the 
baseline and the solar, atmospheric and matter parameters are encoded in these integrals. 
The number of events in the z-th bin, in the same approximation used for eqs. Q, is 
therefore: 



V 



{l+sin 2 (20 13 ) 



iV;+ = {/ 1 -sin 2 (2^ 13 ) + 



1% cos 5 + ig" sin 8 



I 2 cos 5 — J 3 sin 5 



cos^i3 sin (28 13 ) + I A 



(16) 



For a specific energy bin and fixed input parameters 
in the (6*13, 5) plane solving: 



■ 

13, 6), we derive the clone location 



Satm, Soct) 



N l u± (6 13 ,5;s at „ 



±W13 



3, S a tm, S Q ct 1 



S atm, Soct S Q ct) 
^u±(^13; 3] Satm = Satm, S oc t = — Soct) 



N. 



± ($13) S a tm, S ct) 



n: 



±{013 



N, 



,5; s 
,5;s Q 



atm 



±\ u 13, °, Satm 



Satm, Soct 
Satm, Soct 



Soct) 
Soct 1 



(17) 
(18) 
(19) 
(20) 



as it was the case for the transition probability in eqs. ((SUB))- 



3 Location of the clones 

In this section we determine the analytic expression for the clone location, in the case 
of v e — > oscillations, solving the integrated systems in eqs. (jl7H2UJ) . We are going to 
describe completely the procedure for calculating the intrinsic clone and then we present, 
for all the cases, the simplified vacuum formulae. A detailed derivation is deferred to 
Appendix [X] Moreover, in order to simplify the notation, from now on we will omit any 
reference to the bin index % in the convolution integrals, eq. (|15|). 



3.1 The intrinsic clone 

We derive now the theoretical solution of eq. (fTTj) following the outline of [TJ]. This 
equation corresponds to the particular case s a tm = s atm ', s oc t = s oc t- The explicit results 
that we get for 5 and 9 13 can be compared with [7| IT3]. 
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Eq. (|T7|) draws two continuous equal-number-of-events (ENE) lines in the (#13, 8) plane. 
We derive an implicit equation in 8: 



where 



and 



F ± (6) = G ± (9 13 , 0i3, 8) 
F±(8) 



sin 5 



'It 



G ± (9 13 ,9 13 ,8) = 1^1 f(e 13 ,9 1 3)+F ± (8)g(9 13 ,e 13 ) 



(21) 



:22) 



(23) 



The two 0i3-dependent functions, / and g, are: 



f(@i3, #13) 

g(0 13 ,9 13 ) 



sin 2 (20 13 ) -sin 2 (20 



13 J 



cos 0i3 sin(20i 3 ) 
cos 0i3 sin(20i 3 ) 



(24) 



cos 0i3 sin(20 



13 J 



This parametrization is valid for any value of 0i3, including 0i 3 ^ 0: not an irrelevant 
feature, since even in the case of vanishing 0i3 clone solutions with a non null value of 
013 can appear. In this case a totally wrong conclusion about the value of 0i3 can indeed 
originate by a naive analysis of the experimental results. Notice, finally, the trivial limits 
/(013, 0i 3 ) = and #13, 0i 3 ) = 1. 

Solving the system in eq. (j21j) . we derive two equations for sin 5 and cos 8: 



sin<5 = C7(0 13 ,0 13 ) +0(013,013) sm 5 



(25) 



cos 8 

where the coefficients are 

'IU2 



£/(013,013) +0(013, 013) COS 8 



C 



Kit 



UI2 



is4 



D 



I3IJ 



(26) 



In the following results, the baseline and energy bin dependence appears only through 
these integral ratios. Notice that only three of the four integrals if are present in C and 
D. In fact, there is no dependence on I4 due to the fact that for s oc t = s oc t this integral 
drops between left and right hand sides of eq. t)16j). 

In the approximation cos0i3 ~ cos0i3 ~ 1, perfectly valid since 0i3 < 10° jSj, we solve 
eq. (J23J) for i3 to get: 



sin 2 20i3 




D cos 8 + C sin 8 
C 2 + D 2 



sin 20 



13 



(27) 
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in which the first solution represents the true input value #13 = #13 and the second one 
the location of the intrinsic clone. Eventually, substituting the results of eq. (J2*7)l into 
/(#i3, #13) and g(0i3, #13), we can get the value of 5 for the intrinsic clone using one of the 
relations in eq. ()25|). 

It is interesting to specialize these results to the case of very short baseline L, where 
matter effects can be neglected. In the vacuum approximation, the convolution integrals 
of the four Oj coefficients no longer differ for neutrinos and antineutrinos, lim ( 4_ + o if — Ij- 
As a consequence, 

lim C = lim D = — (28) 

and 



sin 5 = #(6>i3, 0i3 ) sin 5 



COS 5 = ^/(013,013) +0(013,013) COS 5 



(29) 



The location of the intrinsic clone in vacuum is: 
sin 2 20 13 = sin 2 20 13 + (j) 



(30) 



Notice that in vacuum the clone location depends on I\ and I2, only: the dependence 
on J 3 has also dropped from eq. (j22J). Eq. (j^Oj) generalizes the results of [7 j and [TB] , 
where two regimes ("atmospheric" and "solar", for large and small 0i3, respectively) are 
discussed independently. In particular, at first order in I2/I1 our vacuum solution for 13 
reproduces eq. (5) of ^3], valid in the atmospheric regime: 

sin 20 13 ~ sin 2 20 13 + 2-^ cos 5 sin20 13 (31) 

h 

cot 5 ~ -coa-^ , 1 . (32) 
li sin 2&i3 sin 

For values of 0i3 large enough, 5 ~ n — 5 + 0{lijl\)- 



3.2 Of other clones 

The general case is represented by eq. (|2*0jl . In this case, integrals on left and right hand 
side of eq. fjlfij) are different and we can no longer factorize easily the 0i3 dependence. 
As a consequence, the notation is going to be more involved with respect to the intrinsic 
case. To gain analytical insight of the results, however, we can carry on a bit more our 
mathematical manipulation. We label as Ij the true integrals (on left hand side of eq. (|16|)) 
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and Ij those to be determined in the experiment (on right hand side of eq. (|16p). We 
modify the functions in eq. (1213) as follows: 



/4(013, 6*13) 

h'(e 13 ) 



sm 2 (26 13 )-[lf/lt} sin 2 (20 13 ) 



(1 - mi 



cos ^13 sin(26>i 3 ) 
1 



(33) 



cos6 l 1 3 sin(26> 13 ) 



with /4 (^i3i ^13) = 7(^13)^13) an d h'(9is) = in the intrinsic case. 
With these definitions we can rewrite eq. (J25J) as follows: 



^(013,013,5) 



I -lme 13 ,9 13 ) + I 'F ± (5) 9 (e 13 ,e 1 3) + ^ti(e 13 ) 



where, clearly, 



F ± (8) = cos5± (b^j sin 5. 



(34) 



(35) 



In Tab. ^ we present the relation between the Ij true integrals and those that undergo 
a discrete s atm and/or s oct transformation. As it can be seen, the effect of s atm —s at m 
is to flip neutrinos with antineutrinos (notice that I4 is unaffected by this change). The 
change in the 9 23 octant introduces a dependence on 623 in Ii and J 4 , leaving J 2 ,3 un- 
modified. Using Tab. ^ we can specialize eq. (JSlj) for the three ambiguities and derive 
analytical results for the sign, octant and mixed clones. These results are typically rather 
cumbersome and we refer the interested reader to Appendix [X] for the explicit analytical 
formulae for the different clone locations. 





Intrinsic: 


Octant: 


Sign: 


Mixed: 




Satm Satm 


Satm S a f m 


Satm S a i m 


Satm S a i m 




Soct S oc f 


Soct S oc i 


Soct S oc f 


Soct S oc i 


It 
It 
It 

n 


It 
It 
It 

n 


cot 2 (9 23 )lt 
It 
It 

tan 2 (0 23 )/ 4 


It 
-It 
It 

n 


cot 2 (0 23 )7f 

-It 
It 

tan 2 (fl 23 ) h 



Table 1: Transformation laws of the integrals L under the discrete ambiguities in the 9 



23 



octant and in the Am 2 , 



23 Sl 9n- 



We specialize here to the vacuum case. In this case, since the convolution integrals 
make no distinction between neutrinos and antineutrinos, the transformation properties 
of Tab. H] simplify considerably. In particular, I 3 is unaffected by any of the discrete 
transformation; J 2 changes sign when s atm —> —s a tm an d is unaffected by an octant shift; 
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I\ and J4 are unaffected by changes of s atm , whereas get a #23 dependence for an s oct 
change. Notice that for #23 = 45°, only J 2 is sensitive to the sign of Am| 3 . 

The vacuum limit of eq. (jHHj) is given by: 



Intrinsic and Sign Clone 

j±(013, ^13) 

< 

. KIM 

Octant and Mixed Clone 

f±(9i 3 ,9i 3 ) — ► fb(9i3, 9i 3 ) -- 



- /(*U,01s) 

-> 

sin 2 (20 13 ) -cot 2 23 sin 2 (2# 13 ) 



(36) 



/l'(^13) 



W13) = (l- tan 2 23 



cos 0i3 sin(20 13 ) 
1 



(37) 



cos 0i3 sin(20i 3 ) 



Eq. (J29J) for sin S becomes in vacuum: 

sin 5 = #(013,013) sin 5 



(3* 



for all the four clones. The equation for cos 5 in vacuum is, on the other hand, different 
in the four cases: 



Intrinsic Clone: 
Sign Clone: 

Octant Clone: 
Mixed Clone: 



(COS S) int = ^ /(0i3, 0i3 ) + #(013, 013 ) COS 5 
(COS S) sign = -(C0s5)int 



(COS S) oct = ^ ^(013, 01 3 ) + #(013, 013) COS 6 + £ h' (9 13 ) 



h 



(cos S) 



mixed 



(cos S) 



oct 



Since the only difference for intrinsic and sign clones and for octant and mixed ones 
is a sign in front of the cos S equation, we only get two vacuum solutions: 



Intrinsic and Sign Clones 



sin 2 20 



13 



sin 2 20i3 



sin 2 20i3 + ( 4- 



2 r 



1 + 2 Ij-j cos6 sin 20i3 



(39) 
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• Octant and Mixed Clones 



sin 2 2# 13 = sin 2 2# 13 + 



(40) 



+ I tan 2 9 23 J- 



± tan 2 9oi -= 



h 



cos 5 sin 2#i3 + tan 823 

2Ii 



— — Q l % 

cos 5 sin 2# 13 + tan # 23 —y- 







tan 2 623) ( sin 2 2#i 3 - tan 2 9 23 j 



- 1 - tan z 



7 23 



sin 2# 13 — tan 9 



2 a J 4 



2:-! 



h 



1/2 



Notice that for the intrinsic and sign clones in vacuum, one solution is the true solution 
#13 = #13 whereas the second solution get shifted with respect to 6>i 3 . This is not the case 
for the octant and mixed clones that, also in the vacuum limit, will undergo a shift with 
respect to #i 3 for both solutions. Eq. can be immediately recovered from eq. (J4*U|) for 
#23 = 45°. 

The solution for the clone location along the 5 axis can be found using eq. ()38|). Notice 
that, due to property of the vacuum transition probabilities under a change in the sign of 
Am 2 3 , the 5 shift for the sign and mixed clones is obtained replacing S sign = 7r — 5 int and 

mixed 7T Ooct- 



3.3 The v e —> v T transition 

The results of Sect. 13.11 and 13.21 can be straightforwardly extended to the v e — > v T oscilla- 
tions. In the same approximation used for eq. (jHJ), we get for the oscillation probability: 

P±(#i 3 , 6) = XI sin 2 (2# 13 ) - Y ± cos(# 13 ) sin(2# 13 ) cos (±5 - ^) + Z T , (41) 
with the following coefficients: 




The coefficient Y± is defined in eq. (II 1J) . and it appears multiplied by a factor (—1) with 
respect to the same term in eq. Q. The different correlation between # 13 and S due to this 
sign has been shown to offer a possible solution of the intrinsic ambiguity at a Neutrino 
Factory, [HI [J3]. 

Eq. (|14j) must be modified to take into account the r — > /i branching ratio and the 
energy distribution of muons proceeding from a r-decay. The number of wrong-sign muons 
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in the detector is therefore: 



dN fi ^(E fl , E T ) d<7 UT {p T )(E T ,E v ) 



dE„ 



dE T 



p±(E u ,e 13 ,5) 



d<t> 



(E„, E^ 



Ei+AE U 



dE u 



The integrals Ij are consequently: 



dN^(E^E T ) da VT(ir) {E T ,E v ) 



dE„ 



dE r 



dE u 



(44) 



The analytic results of Sect. 13. II and 13.21 are formally extended to the v e — > v T oscillation 
replacing the If with IJ as follows: 



If -> tan 2 # 23 /[ 
T 1 * — > -T T 



IS 
i£ 



(45) 



COt 2 023/1- 



4 The Degeneracies flow 



In this section we will apply our analytic results to illustrate, from a theoretical point of 
view, where the different clones of the true solution, (#i 3 ,<5), are located in the (#13, 5) 
plane for different experiments that have been discussed in the literature. The goal 
is to show which combination of experiments is better suited to solve some (or all) of 
the degeneracies. Notice that the existence of unresolved degeneracies could, in fact, 
manifests itself in a complete lost of predictability on the aforementioned parameters. It 
must be noticed that this analysis only provides a useful tool to detect the experiment 
synergies that are most promising to measure these two still unknown entries of the PMNS 
matrix. These results must then be confirmed by a detailed analysis in which statistics 
and systematics of a given experiment combination are carefully taken into account. This 
experimental analysis is well beyond the scope of this paper where we only focus in giving 
a theoretical insight of the problem. 

The experiments that we have considered in our analysis are the following: 

• A Neutrino Factory with 50 GeV muons circulating in the storage ring and two 
detectors: the first one located at L — 2810 Km is designed to look for the golden 
channel (hereafter labeled as NFG); the second one located at L — 732 Km is 
designed to look for the silver channel (hereafter labeled as NFS). The average 
neutrino energy is: < E UejPe >~ 30 GeV. 
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• A SuperBeam with 2 GeV protons and one detector located at L — 130 Km to look 
for — > u e (hereafter labeled as SB). The average neutrino energy is: < E v >~ 270 
MeV and < E H >~ 250 MeV. 

• A Beta Beam with 6 He and 18 Ne ions and one detector located at L — 130 Km 
to look for v e — > (hereafter labeled as BB). The average neutrino energy is: 
< E Ue >~ 360 MeV and < E Dc >~ 240 MeV. 

Notice that for our theoretical analysis the details of the detectors (such as the mass or 
the specific technology adopted to look for a given signal) are irrelevant, since they cancel 
between left and right hand side of eqs. (fTT|) -()20| ) . However, the relevant parameters (the 
neutrino fluxes and the baseline) have been chosen considering the following proposed 
experiments: [1011111 (NFG), [HI EH] (NFS), [U] (SB) and [12 ; (BB). For the different 
beams we have taken as representative the CERN Neutrino Factory, SPL and BetaBeam 
proposals 4 . 

The input parameters used in the computation are: 

• The Atmospheric Sector: Am^ = 2.5 x 10~ 3 eV 2 , 6» 23 = 40°; 

• The Solar Sector: Amj 2 = 7.3 x 10- 5 eV 2 , 6 n = 35°; 

• The Matter Parameter: A = 1.1 x 10" 4 eV 2 /GeV . 

In Sect. 14.11 we compare the clone locations derived from equiprobability and ENE 
curves; in Sect. 14.21 we consider the intrinsic and the sign clones; in Sect. 14.31 we focus 
on the octant and the mixed clones. We have shown in Sect. 13.21 that the two pairs of 
solutions are indeed related in the vacuum limit. 

4.1 Equiprobability flow vs ENE flow 

From a theoretical point of view one can look at the clone displacement with respect to 
the true solution (#i 3 = 0i 3 , 8 — S) solving either the systems for the equiprobability 
curves, eqs. ([HHBJ), or the corresponding systems for the ENE curves, eqs. (|17H20|) . One 
should expect the two treatments to be equivalent: this is not always the case. 

In Fig. n we plot the location of the intrinsic clone that results from the equiprobability 
curves and from the ENE curves. These points graphically depict how far from the true 
solution the degeneracy falls. A small change in the input parameter #i 3 results in a small 
shift of the clone location (remind, however, the 27r-periodicity in the 5 axis). Almost 
continuous geometrical regions where degeneracies lie are therefore defined for a given 
interval in # 13 , illustrating how the clones move due to a change in the input parameters: 
we will call this the "clone flow". In the left plot we present the results for the Neutrino 

4 See for example in [21] for a detailed description of each of these proposals. 
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Figure 1: Intrinsic clone flow for the Neutrino Factory golden and silver channels 
(left) and the SuperBeam and BetaBeam signals (right), in the (A9i 3 ,5) plane for 
9 is G [0.1°, 10°], 5 = 90°. The full circle is the true solution. The full (empty) sym- 
bols represent the intrinsic clone flow, computed from the ENE (equiprobability) curves. 



Factory-based golden and silver channels (i.e. v e — > u^, u T ), whereas in the right plot those 
of a SuperBeam-based — > v e and of a BetaBeam-based v e — > signals. The clone 
locations have been computed for variable input parameter 9\% in the range #13 G [0.1°, 10°] 
and fixed 5 = 90°. The arrows describe the direction of the flows: from large to small ^13. 
The plots should be read as follows: 

• the horizontal axis is A#i3 = 9i3 — #13. This is the value of the ^-shift, as derived 
in the second solution of eq. (|27jl : 

• the vertical axis is the CP-violating phase S itself, as it can be derived introducing 
the result of eq. (|27j) into one of eqs. 

For any value of the input parameter #13 the true solution is represented by the full 
(black) circle in A9i% = 0, 5 = 5 = 90°. Notice that this property is specific to the 
intrinsic case. In general, both solutions of eqs. (jl8j) - PUj) will differ from the true solution 
(#13,(5). The intrinsic clones, computed using the ENE or the equiprobability equations, 
are respectively plotted as full or empty symbols as explained in the legend and in the 
caption. 

In the left plot of Fig. can clearly identify the golden and silver channel intrinsic 
clone flow 5 . The NFG flow goes from the true solution toward positive A9\3 and positive 

5 For this specific plot we are choosing the energy bin between 30-40 GeV but the same conclusion 
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5-5: for #13 = 10° we have A9 13 ~ 0.5° and for very small 9 13 , A9 13 -> 1.5°. The 
NFS flow goes from negative A# 13 and positive 5 — 5 toward positive A# 13 and negative 
5 — 5. The different behavior of the two flows reflects the results of eq. (|27|). For u e — ► z^ t 
transitions, the second solution of eq. (J27)l gives sin 2 2# 13 > sin 2 2# 13 for any value of 9 i3 , 
being C and D strictly positive quantities. On the contrary, due to the formal replacement 
of the convolution integrals 1% — > — JJ,/ 3 — > —I 3 we find that C T , D T are strictly negative 
quantities. For changing 9\ 3) therefore, the #i 3 -shift can either be positive or negative. 
Notice how the intrinsic flows computed from the integrated ENE curves and from the 
equiprobability curves (almost) coincide for the Neutrino Factory-based samples. 

In the right plot of Fig. SuperBeam and BetaBeam-based flows are shown. The 
ENE clone flow for the SB facility (full diamonds) has a very small A6* 13 (< 0.1°) and it 
is practically insensitive to the change of 9\ 3 . The shift in 5 — 5 is always negative and 
become smaller for larger 9\ 3 . The ENE clone flow for the BB facility (full squares) is 
practically identical to the SB one. This means that the presence of the intrinsic clone 
in these two facilities does not interfere with the measure of 9i 3 , whereas can drastically 
reduce the ability of measuring the CP phase. 

The important feature is that for all these facilities the difference between using ENE 
(full symbols) or equiprobability (empty symbols) curves is quite substantial. Using the- 
oretical information from the equiprobability curves could lead to largely wrong conjec- 
tures, in particular regarding the 9\ 3 dependence that is practically absent in the ENE case 
and is conversely quite large using the equiprobability curves. The reason why SB and 
BB present this large difference between the two treatments has to be searched in their 
quite narrow energy spectrum. For very peaked spectra the approximation of "constant 
flux" , that is implicit in the equiprobability curves, badly fails. On the other hand, the 
broad neutrino spectrum from a Neutrino Factory makes the two treatments practically 
equivalent. 

The theoretical results obtained for the intrinsic clone are quite general: a large differ- 
ence in computing the clone location from ENE or equiprobability curves is observed for 
narrow band beam experiments, whereas no difference is observed for broad band beam 
experiments. For this reason, in the following sections we will compute the clone location 
for the different experiments using the integrated ENE curves. 

4.2 The Intrinsic and the Sign Clone Flow 

We now look for the different combinations of proposed facilities that help in solving 
the eightfold-degeneracy. The theoretical requirement for solving ambiguities is to find 
a specific combination of experiments such that the corresponding clone flows lie "well 
apart" . In this scenario an experimental fit to the data will result in a good x 2 absolute 
minimum near the true solution, where the information from different experiments adds 
coherently. Near the theoretical clone locations poor x 2 relative minima will be found, 

holds for all the other bins as well for the integrated flux. 



14 



■ BB 

♦ SB 

* NFG 
a NFS 

• TRUE 



-0.5 



0.5 
Afln 



1.5 



180 



135 



-45 



-90 



-135 



If/ 



■ BB 

« SB 

* NFG 
a NFS 

• TRUE 

-0.5 0.5 

Aft, 



1.5 











* 






• 










♦ 


A 








H 






■ 


BB 








♦ 


B+SB 








★ 


NFG 








A 


NFS 








• 


TRUE 









-0.5 



0.5 
Aft, 



1.5 



180 



135 



<5 CP 



-45 



-90 



- 1 35 















•J 






A 






■ 


BB 






« 


SB 






* 


NFG 






A 


NFS 






• 


TRUE 












A 



-0.5 



0.5 

Aft 3 



1.5 



Figure 2: Clone location in the (A9is,5) plane for the intrinsic and sign degeneracies: 
BB (boxes); SB (diamonds); NFG (stars) and NFS (triangles). The full circle is the true 
solution, A6 13 = 0, 5 = 90°. Pictures (a) and (b) represent the intrinsic and sign clone 
flows for 6 13 G [0.1°, 10°] and fixed 5. Pictures (c) and (d) represent the intrinsic and sign 
clone location for 6 13 = 1° . 



since in these points the results of different experiments do not add coherently. In Fig.|21we 
present the clone flow for the intrinsic (left) and the sign (right) degeneracies. The clone 
locations have been computed for variable input parameter ^13 in the range #13 G [0.1°, 10°] 
and fixed 5 = 90° (upper plots) or for a fixed value, #13 = 1°, 8 = 90° (lower plots). Again 
the arrows indicate the direction of the flows from large to small ^13. In each figure the 
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results for the four facilities are presented together, in order to show how the combination 
of two or more of them can help solving that specific degeneracy. We illustrate the outcome 
of our analysis by studying the two degenerations one by one. 

• The Intrinsic Clone: Fig. Et/Efc- 

The SB and BB intrinsic clone flows are extremely similar. This has already been 
noticed in the previous section, and is the result of the beam design. In [TS] has 
been proposed to run the SPL and the BetaBeam at the same time and using the 
same detector located at L = 130 Km. In view of solving the degeneracies, this is 
clearly not a good option since the two flows always lie very near. The combination 
of these two facilities will still be affected by the intrinsic ambiguity. 

This is not the case for the Neutrino Factory setup: the NFG channel is well sep- 
arated from the NFS channel for any value of #13. The golden channel combined 
with the silver and/or the SB/BB facilities is always, in principle, capable to solve 
the intrinsic ambiguity (i.e. provided that the corresponding experimental samples 
are all statistically significant and that systematic errors under control). 

In Fig. |2t we plot the intrinsic clone location for the small value # 13 = 1°, on the 
verge of the sensitivity limit. 

• The Sign Clone: Fig. ^/Hi- 
It can be seen that again the SB and BB flows lie very near. A combination of 
them is therefore not suited to solve the sign degeneracy. On the other hand the 
Golden and Silver Neutrino Factory flows are quite separated. There could be some 
superposition between all the four flows for certain values of #13 (e.g. NFG and 
NFS flows in the limit of vanishing # 13 ). Nevertheless this superposition occurs in a 
region nearby the true solution where the existence of a clone is somewhat irrelevant. 
The combination of either NFG or NFS with the SB or the BB signals should easily 
solve the ambiguity 6 . In Fig. |2t we plot the sign clone location for #13 = 1°. Notice 
that for this specific case no NFG clones were found. This can be interpreted as 
a consequence of the very good (theoretical) sensitivity to matter effects (i.e. the 
Am| 3 sign) of the NFG for L = 2810 Km. Be aware, however, that in a realistic 
experimental case a sign clone could appear also for the NFG due to statistical and 
systematical uncertainties. 

As a final comment notice that the results for the sign clone for the SB and BB flows 
are in perfect agreement with the theoretical analysis of Sect. 13.21 For these two facilities 
the vacuum limit can be considered a good approximation (being L = 130 Km) and we 
therefore expect the same A9i 3 shift as in the intrinsic case with <5 sign = n — 5 int . This is 
indeed the caseby looking at Fig. |2t and Fig. 

6 Notice, however, that for large #13 angles even the Golden channel alone can solve the sign ambiguity 
being strongly sensitive to matter effects. 
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Figure 3: Clone location in the (A6i 3 ,5) plane for the octant and mixed degeneracies: 
BB (boxes); SB (diamonds); NFG (stars) and NFS (triangles). The full circle is the true 
solution, A9i 3 — 0,5 — 90°. Pictures (a) and (b) represent the octant and mixed clone 
flows for 6 13 G [0.1°, 10°] and fixed 5. Pictures (c) and (d) represent the octant and mixed 
clone location for 9i 3 = 1°. The results have been computed with 6 23 = 40°. 
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4.3 The Octant and the Mixed Clones 



If 023 was maximal no additional degeneracies were present. Otherwise two furtherde- 
generacies do appear 7 . In Fig. El we present the clone flow for the octant (left) and the 
mixed (right) degeneracies. The clone locations have been computed for variable input 
parameter 13 in the range #13 G [0.1°, 10°] and fixed 5 = 90° (upper plot) or for a fixed 
value, ^13 = 1°,5 = 90° (lower plots). In each figure the results for the four facilities are 
presented together. 

In Figs. Ek and Eb a very complicated pattern of clone flows appears. For this reason 
we discuss the specific case of 6*13 = 1° depicted in Figs. Efc and Eli. 

• The Mixed Clone: Fig.Eb/Eli. 

We again notice that BB (boxes) and SB (diamonds) clones lie very near. As in the 
case of the sign ambiguity, no NFG mixed clones appear for this specific #13 . Finally, 
the NFS clones are well separated from the BB/SB pairs. The combination of the 
golden and/or the silver channel and one of the BB/SB facilities can, in principle, 
solve the mixed degeneracy. 

• The Octant Clone: Fig. Ek/Efc- 

This is certainly the most difficult degeneracy to solve, since the flows of the four 
facilities overlap in a significant region of the (A# 13 ,(5) plane, Fig. EK- It seems 
feasible to reduce the overall ambiguity by combining NFG and NFS or BB/SB 
facilities, as it can be seen in Fig. Eh- A complete elimination of this degeneracy can 
perhaps be achieved combining the NFG, NFS and one of the BB/SB experiments. 

5 Conclusions 

In this paper we derive the theoretical location of the so-called "clones" that appear 
in the simultaneous measurement of #13 and 5. The general technique used to write 
analytical formulae for the clone location for the four degeneracies (called intrinsic, sign, 
octant and mixed) has been reported. Our analytical results apply to equiprobability and 
equal-number-of-events (ENE) curves as well. However, we have shown that the clone 
location computed from equiprobability or ENE curves can differ significantly for narrow 
band beam experiments. Notice that the latter approach reproduces the experimental 
situation, where the signal is represented by charged leptons in the detector and not by 
the oscillation probability itself. Using theoretical information from the equiprobability 
curves could thus lead to largely wrong conjectures. In the actual calculation we have 
therefore presented results obtained using the ENE curves. 

7 Be aware that these clones will appear also for maximal #23 once including the unavoidable experi- 
mental error band. 
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We present the clone flow for experiments fuelled by three different beams, using as 
a reference the CERN proposal for a SuperBeam, a BetaBeam and a Neutrino Factory. 
In the case of the Neutrino Factory two different detectors have been considered to take 
care of the "golden" and "silver" channel. 

Our results show that the combination of these specific SuperBeam and BetaBeam 
does not help in solving any of the degeneracies. All the clone solutions for these two 
facilities lie very near, as a consequence of the fact that they have a very similar energy 
spectrum and that both are using the same detector located at L — 130 Km. 

The combination of the Neutrino Factory golden and silver channel, on the contrary, 
seems to help significantly to reduce ambiguities. In our theoretical analysis, indeed, the 
intrinsic, sign and mixed clone flows for NFG and NFS lie well apart for d\ 3 > 1°. As a 
general comment, we deduce that the combination of the NF golden channel with either 
the NF silver channel or with one of the pair BB/SB could in principle solve these three 
ambiguities. 

This is not the case for the octant clones. Although a significant reduction of the 
degeneracy can be attained combining the NFG with the NFS or with one of BB/SB, it 
seems that a complete elimination of the octant clones can be achieved only by combining 
the three facilities, NFG, NFS and BB/SB. Nevertheless, even using the combination of 
all the available facilities is quite possible that the experimental uncertainties will not 
permit to solve completely the octant ambiguity for small values of #13. The ultimate 
sensitivity of the future neutrino experiments to S is therefore strongly dependent on the 
success or failure in solving it. To decouple the pair BB/SB (e.g. using the Minos OA 
proposal (221 f° r a SuperBeam) can probably be of great help in solving the octant clone. 
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A The Analytic Clone Location in Matter 



In this appendix we generalize the treatment given for the intrinsic clone in Sect. 13. II to a 
generic case where the discrete variables s atm and s oct can assume different values between 
left and right hand side of eqs. (JT7j)-(j2nj). 

Using eq. (|3^j) and eq. (|3~2j) we get: 



sin 5 



C\ sin 26 13 - d sin 2 26 13 



cos 8i 3 sin 26 
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Co cos 5 
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7 13j 
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(46) 
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D\ sin 29 13 - D 1 sin 2 20i3 
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D' 2 cos 5 + C 3 sin 6 g(9 13 , 9 13 ) + D' 4 h'(9 13 ) 



cos 9 is sin 2#i3 

where, using a notation as much resembling as possible that of the intrinsic case, we have 
defined the following coefficients: 
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Notice that in vacuum all the coefficients of the C-type vanish. Using the discrete 
transformation properties of the integrals Ij of Tab. [I]we derive the corresponding trans- 
formation properties for the Cj and Di coefficients, reported in Tab. El The 2 3 -octant 
transformation does not affect any of the C[ or D[ coefficients, but only C\ and Di that are 
multiplied by cot 2 6*23 when 6*23 — > n/2 — 9 23 . The effect of an s atm change is much stronger: 
all coefficients are affected by this change. The C[ and D[ coefficients are obtained by the 
Ci, 2 ,3> £>i,2,3 in eq- (S3) b y changing if -> -1$ and if -> 7 3 T . 

We present hereafter the explicit analytic solution for the #i 3 -shift for the octant and 
sign clones. Inserting these results in one of eqs. (J4lj)) the 5-shift can be extracted. We do 
not present here the corresponding expression for the mixed clone. 



A.l The Octant Clone 

Using Tab. |2l we get for the #i 3 -shift for the octant clone the following expression: 
sin 2 2fl 13 = sin 2 2# 13 + (4f 
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Table 2: Transformation laws of the coefficients Ci,C[ and Di, D[ under the discrete 
ambiguities in the #23 octant and in the Am 2 3 sign. 
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It can be easily shown that in the vacuum limit, when C\^C[ — > 0, this equation 
reduces to the vacuum solution, eq. (jlOJ). Moreover, for 2 3 — > 45° reduces to the general 
solution in matter, eq. ()27j) . for the intrinsic clone. 



A.2 The Sign Clone 

Using Tab. [21 we get for the 0i3-shift for the sign clone the following expression: 
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((7iC£ + D^) 2 + ((7^ + D 1 D[ 
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-8 (CiD; - DiC^) [(Ci^ - DiCj) cos 5+ (C1C3 - £>i£> 3 ) sin 8] sin 3 2fl 13 
-4 (pxD' 2 + Dx^) 2 cos 2 5 sin 2 2# 13 
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One can verify that eq. (jlQj) reduces to the vacuum solution, eq. (JSHJ- The first term 
in the square root goes to a perfect square; the second vanishes for an exact cancellation 
between (DiD' 3 ) 2 and DiD[, since in vacuum D[ — > — Z)^ — > —1 and Z) 3 — > 1. All other 
terms in the square root trivially vanish since are multiplied by some C-type coefficient 
that goes to zero in the vacuum limit. 
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